Origin of the anomalous Slater-Pauling curve in CoMn alloy 

clusters 

M. Pereiro/'Q D. Baldomir/ and J. E. Arias^ 

^ Departamento de Fisica Aplicada and Instituto de Investigacions Tecnoldxicas, 
Universidade de Santiago de Compostela, 
Santiago de Compostela E- 15782, Spain 
^Instituto de Investigacions Tecnoldxicas, 
Universidade de Santiago de Compostela, 
Santiago de Compostela E- 15782, Spain 
(Dated: August 17, 2009) 

Abstract 

Surprising enhancement of the magnetic moments recently observed in dilute CoMn alloy clusters 
is explained using ah initio electronic structure calculations. The calculated magnetic moments 
generally agree with the reported experimental data. An equation for calculating the magnetic 
moments of the CoMn alloy clusters has been derived to correct the deviations predicted by the 
rigid-band model and the virtual bound states approximation. A new strategy is proposed to 
obtain the ground-state structures of the CoMn clusters and it was also put to the test of the 
experiment. 

PACS numbers: 36.40.Cg,75.20.Hr,31.15.A- 
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I. INTRODUCTION 



The exploration of bimetallic transition metal (TM) clusters is emerging as a promise 
field of research because of the new opportunities they offer for developing magnetic record- 
ing devices and cluster-assembled materials with new functions for medical applicationsii^. 
Although the electronic and magnetic properties of bare TM clusters have been actively 
studied for several years, however less attention has been paid to their alloys because they 
represent both an experimental and theoretical challenge. Thus, for experimentalists is 
very difficult to control the stoichiometry of the alloy clusters using chemical methods, and 
for theoreticians, the determination of the ground-state geometries becomes a very difficult 
task, as commented later in the main text. The recent observation of the average magnetic 
moment enhancement in dilute CoMn alloy clusters^ that contradicts the bulk behavior has 
introduced an entirely unexpected dimension on the subfield of magnetic alloy clusters and 
it could pave the way for future possible applications like, for example, in biomedicine as 
magnetic sensors. Although the search for the origin of the aforementioned observation has 
actively stimulated the research on CoMn alloy clusters^i^i^i^, however the answer still re- 
mains elusive. A possible explanation is based on the assumption that a virtual bound state 
(VBS) is formed below the highest occupied molecular orbital (HOMO) level and near the 
Mn site^. A VBS can be defined in the potential scattering model as a resonant scatter- 
ing near impurity atoms in the host which induces a narrow peak in the conduction band 
density of states (DOS). Originally, the VBS model was developed by Friedef'' to explain 
many of the physical properties of bulk magnetic alloys containing dilute magnetic species. 
It represents an improvement over the rigid-band (RB) model which is based on the assump- 
tion that the s and d bands are rigid in shape as atomic number of the alloy changes^. In 
this article, our calculations demonstrate that the origin of the anomalous behavior of the 
Slater-Pauling (SP) curve of CoMn alloy clusters does not require the formation of a VBS 
as suggested in Ref. Isl, but is explained directly in terms of the magnetic moment provided 
by the Mn atoms and the "spin-flipping" of the electrons belonging to the Co-Mn bond- 
ing. Thus, with the aim of studying the magnetic properties of Co-Mn alloy clusters and 
elucidating the anomalous behavior of their SP curve, we have performed spin-unrestricted 
density functional theory-based (DFT) calculations as implemented in the demon-ks3p5 
program package^. 
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II. METHOD AND STRUCTURAL PROPERTIES 



The electronic system consisting in nuclei and r] electrons is assumed to be described by 
the next Hamiltonian, which expressed in second quantization and in the Born-Oppenheimer 
nonrelativistic approximation reads as 

/ 1 ^ \ 

ij ffi \ k ) 

1 

ijkl ufiaf3 

(p^ ^pn (1) 

where the Kohn-Sham orbitals were expanded into atomic orbitals \i) = Y^^c^i in 
the Linear Combination of Gaussian Type-Orbitals ansatz. The orbital basis sets of con- 
traction pattern (2111/211*/311+) and (2211/311/411) were used in conjunction with the 
corresponding (5,3,3;5,3,3) and (3,4;3,4) auxihary basis sets for Co and Mn, respectively^. 
Likewise, ad hoc model core potentials with contraction pattern (4:7,4) and (5:7,4) have 
been used for describing the inner electrons of Co and Mn atoms, respectively^. The ki- 
netic and nuclear attraction energy of the electrons in an environment of Zk nuclear charges 
are described by the one-electron term. The || symbol represents the l/|r — k| operator, 
where r and k are the electron and nuclear position vectors, respectively. The two-electron 
operator represents the Coulomb repulsion energy of the electrons. In this case, the symbol 
II represents the l/|ri — r2| operator for the electrons with coordinates ri and T2- The last 
term in Eq. ([1]) is the exchange-correlation energy and we have used here the form proposed 



in Ref . 



12 which is a function of the spin-dependent electron density (p'^, cr =],[) and its gra- 



dient in the generalized gradient approximation (GGA). We have adop 
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;ed the former GGA 
that it represents a 



functional in our DFT calculations because it was reported in Ref. 
dramatical improvement in the calculated binding energies of Co clusters. The total energy 
of the system is calculated adding the nuclear repulsion energy {J^^^bi^aZb) /\^a — '*^b\) to the 
electronic contribution. A wide set of spin multiplicities ranging from 1 up to a maximum of 
61, depending on the selected cluster, was checked to ensure that the lowest-energy electronic 
and magnetic configuration is reached. More information about the computational details 
can be found elsewhere^. Hereafter, all calculated results refer only to the Co-Mn alloy 
20 atom clusters. We have also performed electronic structure calculations of some guessed 
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geometries of the Co-Mn alloy clusters with Q = 25 and 30 in a range of Mn concentration 
less than 0.3, but the obtained results slightly differ from that of the Q = 20 case, i.e., we 
observed that the calculated magnetic moments of clusters with Q = 25 and 30 increase 
with the impurity concentration as i7 = 20 clusters do. Accordingly, the conclusions drawn 
from the Co-Mn clusters with Q = 20 can be extended to clusters with greater size. 

The search for the global minima of the Co-Mn alloys 20 atom clusters was planned as a 
multistage strategy combining an unbiased search method i.e., a basin-hopping^^ algorithm, 
in conjuction with a molecular mechanics methoci^. In the reoptimization procedure of the 
clusters, we have made use of the Polak-Ribiere algorithm^ without any symmetry con- 
straint and the root-mean-square gradient was set to 10^'^ Kcal/(mol A). In a first stage of 
the calculation method, the initial guessed structures of the C020 cluster were taken from 
three different sources, namely, the structures were provided by the GMIN code'^^ which uses 
the basin-hoping algorithm, the existing databasesi^, and proposed by us. After that, they 
were reoptimized with the HyperChem codecs with the intention to obtain the ideal candi- 
date to the lowest-energy structure of the C020 cluster. The converged guessed structures 
are illustrated in Fig. [1] and the structural and energetic parameters are reported in Table [B 
The ground-state structure (Co2o_o) is a capped double icosahedron. The cohesive energy 
and the number of nearest-neighbor Co-Co bonds reveal that the most stable structure is 
also the most compact one although the average first-neighbor distance is a little bit higher 
than the rest of structures plotted in Fig. [H In a second stage, the rest of guessed structures 
of the ComMn„ alloy clusters were built from the geometry of the C020 cluster but this time 
m — n Co atoms swap positions with n Mn atoms. The assumption is underpinned by three 
arguments: firstly, the ionic radii (r^) of Co and Mn atoms are pretty close i. e., 1.25 A and 
1.26 A, respectively^. Secondly, we have calculated the size of the Thomas-Fermi screening^ 
outside the atomic sphere of the Mn atoms (Zmh) 

Ztf = Zmu (1 + ra^^nr]{eHOMo)) c^'-'^V^-'^I^homo) (2) 

and it is vanishingly small because the DOS at the HOMO level {rj^enoMo)) provided 
by our DFT calculations are large [see Fig. [2](c)], so that the exponent in Eq. ([2]) is 
ra^j4:Tr {rj{eHOMo)) — 247.58, where {rj^enoMo)) represents the average over the lowest- 
energy Co-Mn alloy series clusters with Q = 20. Thus, the size effects are negligible and 
the addition of Mn atoms to the cluster should slightly distort their environment. Thirdly, 
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FIG. 1: (Color online) Illustration of the guessed structures of the C020 cluster. The structures are 
ordered from left to right and top to bottom by increased relative energy. The notation Co2o_m 
stands for the mth energetic isomer. 



we have performed a geometry optimization of C03, Co2Mni, C04, and CosMni clusters by 
means of spin-unrestricted DFT calculations^ and the results show that the elongation of 
the Co-Mn average bond length with respect to the Co-Co one is less than 0.15 A. In conse- 
quence, the substitution of cobalt atoms by Mn atoms should slightly distort the geometric 
structure of the clusters. This fact has also been observed in Ref. sj. 

Another specific difficulty of alloy clusters is the great number of possibilities they offer 
to distribute the solute atoms over the sites at a given geometrical arrangement of the 
atoms. Thus, for example, in the case of CoioMnio the number of possible homotops for a 
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TABLE I: Average first-neighbor distance (dco-Co)i number of nearest-neiglibor Co-Co bonds 
(A'^Co-Co) per atom, relative (AS), and coliesive (Ec) energies of C020 cluster isomers. The geom- 



etry notation is that of Fig. [TJ The point groups are determined from Ref. 



1^. 
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C020_5 
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2.44 


1.50 
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1.38 



geometry with inequivalent positions is C^q = ~ 1.8 x 10^. With the aim of overcoming 
the aforementioned point, we have planned out a new strategy consisting in optimize the 
geometries of the 20 homotops belonging to the CoigMni cluster. As starting point, every 
structure was built from the optimized geometry of C020 cluster but substituting a Co by a 
Mn atom at the 20 accessible positions of the cluster. After that, we performed a geometry 
optimization of the 20 homotops to obtain the stability of them as a function of the position 
occupied by the Mn atom. The rest of Co-Mn alloy clusters were created by replacing the 
cobalt atoms by Mn atoms at the positions with the lowest energy and once again they were 
reoptimized in geometry. We examined the method for the case of Coi8Mn2 homotops where 
one of the two Mn atoms was attached to the position with the lowest-energy and the other 
Mn was positioned in the rest of available sites. After the optimization procedure, we have 
observed that the lowest-energy structure was the one with Mn atoms in the lowest-energy 
positions and the energy of the other geometries was estimated to be 3.22 eV greater than 
the ground-state cluster. We have also put the geometry optimization method to the test 
of the experimenti^S as shown in Table [TTl The ionization potentials provided by our DFT 
calculations are in good agreement with the available experimental results and consequently, 
our reported structures should be very close to the ground-state structures of the Co-Mn 
alloy series clusters with Q = 20. The lowest-energy structures are illustrated in Fig. [3] and 
the geometrical properties are collected in Table IIIII The Co-Mn alloy geometries converged 
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in slightly distorted structures of that of the Co2o- In particular, the distance of the atom 
labeled as 12 in Fig. [3] to the rest of atoms belonging to the cluster is gradually reduced 
with the enhancement of the Mn concentration up to the CoisMus cluster which elongates 
the distance. For the rest of clusters, the distance persists approximately unaltered. This 
point is reflected in the average nearest neighbor distance for Co-Co bonding reported in 
Table IIIII It is also important to comment that the average nearest neighbor distances for 
the Co-Co, Mn-Mn and Co-Mn bondings are in general very close each other (see Table [In|) . 
This results favor the assumption commented above that the addition of Mn atoms just only 
alters a little bit the geometry of the ground-state structure of C020 cluster. To conclude 
the structural discussion, and with the aim to show how the Mn atoms diffuse into the Co 
ones, it is convenient to define an order parameter that is positive and close to 1 when the 
phase separation or segregation takes place and close to when the mixing or disorder is 
the main contribution to the arrangement of the Mn atoms in the cluster. The parameter 
which meets the aforementioned conditions is the chemical order and it is thus defined as 

where Ni_j represents the number of nearest neighbor i-j bonds (with i=Co and j=Mn). 
The values of F reported in Table UTTl are plotted in Fig. HI It is easily observed in Fig. H] the 
asymmetry of the chemical order parameter with respect the dash vertical line. Thus, at 
low concentrations of Mn atoms F decreases more rapidly than for higher concentrations of 
them so that the mixing is favored in an interval of Xmh ranging approximately from 0.15 to 
0.65. The reason for that behavior is attributed to the positions occupied by the Mn atoms 
in the range of low concentration. As can be seen from Fig. [HI the impurity atoms swap 
positions with the inner Co atoms. The number of bondings between the inner positions 
and the rest of atoms is higher than the outer-shell positions. Thus, the number of Co-Mn 
bondings is favored for lower concentrations of the impurity and in consequence F decrease 
rapidly. But however, for higher concentrations of the Mn atoms, the process is not inverted 
and the cobalt atoms occupy mainly the outer-shell positions instead of the inner ones. 
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TABLE II: Ionization potentials (in eV) for ComMn„ clusters with O = 20. The experimental data 



were taken from Ref. 
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5.37 


5.34 


5.31 
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III. RESULTS AND DISCUSSION 

The anomalous behavior of the SP curve in Co-Mn alloy clusters compared to the CoMn 
bulk— is shown in Fig. [51 The magnetic moments per atom (empty star symbol) increase 
with a slope of 1.83 /i^ up to a Mn concentration of 40%, after which the average magnetic 
moments tend to decrease with increasing Mn concentration. The above results are in very 
good agreement with the available experimental measurements (filled star symbol) reported 
in Ref. y (1.7 and 40%, respectively). The discrepancy with the numerical data (about 
0.6 /i_B in average) may be related to the omission of the orbital moment contribution since 
the total magnetic moment is (M) = 2(5*) + (L). However, according to the results of the 
calculations performed for binary TM clusters reported in Ref. Ij, the orbital moments rep- 
resent a very small correction to the total magnetic moment and consequently underpins our 
approximation. Likewise, the source of error could also be ascribed to a wrong assignment 
of the ground-state structures of Co-Mn alloy cluster, but however it is unlikely because the 
measured value of the magnetic moment^^ in the case of the C020 cluster is about 2.04 /i^ 
which is very close to our predicted value of 2.00 hb- 

The negative slope (~ 6.0 /i^ per Mn substitution) of the SP curve for the bulk CoMn 
alloy, which is in contrast to the cluster behavior, is explained on the basis of the VBS 
concept''', i. e., a resonant scattering of the conduction electrons at the Mn sites which 
induces a narrow peak above the Fermi level for both spin channels. This point was confirmed 
by ah initio band structure calculations^*^. In the case of the CoMn alloy clusters, the first 
attempt to explain the positive slope was also based on the assumption of the existence of a 
VBS but this time only the majority-spin VBS was conjectured to remain below the HOMO 
levels'. However, we can rule out the existence of the VBS and come to the same conclusion 
if we admit two hypothesis: first, the assumptions of the naive RB model are valid and 
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TABLE III: Structural properties of the lowest-energy Co-Mn alloy clusters with = 20. We 
report the average nearest neighbor distance for the Co-Co, Mn-Mn and Co-Mn bondings. The 
chemical order defined in Eq. ([3]) is also provided in the last column. The interatomic distances 
are given in A. 
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where c is the impurity concentration per atom and IS.Z is the atomic number difference of 
the impurity relative to the host. For Co host doped with Mn impurities, AZ = —2 and the 
slope of the average alloy magnetic moment relative to the host is positive and proportional 
to 2 yU^, which is relatively close to the experimental value (1.7 /i^) and our reported result 
(1.83 /is). The VBS approximation can predict successfully the SP curve behavior in most 
of the TM solid alloys whereas its application to alloy clusters is less reliable because of the 
lack of a periodic crystalline lattice^^. Moreover, the VBS and RB model fail to predict the 
behavior of the SP curve for Mn concentrations greater than 40 % because they are only 
valid for small concentrations of impurities. 

We have derived a formula for the average magnetic moment of CoMn alloys that corrects 
Eq. (jlj) and takes into account the antiferromagnetic (AF) coupling between the impurities: 

(/^aiioy) = m {fico) + {n- g{n)) - g{n) (/ij^^) (5) 

where m and n are the number of Co and Mn atoms in the cluster, respectively and g{n) 
is the number of Mn atoms that couple antiferromagnetically with the rest of atoms^i (see 
Table IIV|) . In Fig. [3l we have also plotted in green the Mn atoms that couple antiferromag- 
netically to their neighbors. The averaged magnetic moments per atom and per clusters 
with different Mn compositions obtained from our DFT calculations for Co and Mn atoms 
are (/ico) = 1.93, <^/iMn^ = 3.30, and <^/iMn) = 2.15 HB/atom. The values (diamond points) 
provided by Eq. ([5]) are plotted in Fig. [5l To gain more insight into the physics behind 
Eq. ([5]) and based on the results provided by the electronic structure calculations, we have 
identified a double mechanism that explains the magnetic enhancement trend for CoMn 
clusters below 40% of Mn concentration and the successive fast dropping of the magnetic 
moment. 

The first mechanism resides in the role played by the Mn atoms, as reflected in Eq. (JSj). 
First, the clusters increase their magnetic moments due to the addition of the Mn moments 
up to a maximum of 2.7 fiB/O'tom and then they decrease their moments because of the 
AF alignment of some Mn atoms with the rest of atoms belonging to the cluster. We have 
plotted in Fig. [6]the magnetization density of Coi2Mn8 and ConMug to see how the clusters 
evolve from a ferromagnetic (FM) configuration to an AF at a critical concentration of the 
impurity. The red surface surrounding the Mn atom indicates an AF coupling with the rest 
of atoms while the blue surface indicates a FM alignment. The onset of Mn atoms with 
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TABLE IV: Average magnetic moments (in fiB/o-tom) and number of Mn atoms (g(n)) that couple 
antiferromagnetically with the rest of atoms belonging to each cluster of the ComMn„ series, and 
with n ranging from 9 up to 20. The notation N}^^_j^j^ (^Mn-Mn) represents the average number 
of nearest neighbor bonds between a FM (AF) Mn atom and the rest of Mn atoms. 



n 


9 


10 


11 


12 


13 


14 


15 


16 


17 


18 


19 


20 


g(n) 


-0.98 
1 

4.6 
7.0 


-1.11 
1 

4.9 
8.0 


-1.45 
1 

4.4 
8.0 


-1.50 
1 

5.1 
8.0 


-2.39 
3 

4.9 
8.3 


-2.49 
3 

5.3 
8.7 


-2.31 

4 

5.1 

8.5 


-2.73 

5 

5.3 
8.4 


-2.67 

5 

5.7 
8.8 


-2.27 

5 

5.9 
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6.8 
7.3 


-2.95 

8 

6.2 
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negative magnetic moment resides mainly in a Mn-Mn charge transferring. The values of 
^Mn-Mn ^Mn-Mn reported in Table HVl indicate that the AF Mn atoms are the ones which 
establish more bondings with the rest of Mn atoms in relation to the FM Mn atoms. To 
exemplify the explanation, we have plotted in Fig. [7] the spin-polarized DOS for the Mn-13 
and Mn-16 atoms of the ConMng cluster. The Mn-13 atom is coupled antiferromagnetically 
to the rest of Mn and Co atoms (see Fig. [3]). The elevated number of bondings make that 
the Mn-13 atom share more electrons with its Mn environment. This produces a charge 
transferring from the spin- up channel of the Mn-13 atom to the spin- up channel of the Mn 
atoms surrounding it, that is for example the case of Mn-16. Thus, at a critical number of 
the Mn-Mn bondings (close to 7) the Mn-13 atom becomes antiferromagnetic. 

The second mechanism involves also the Co atoms and the "spin-flipping" of the electrons 
belonging to the Co-Mn bonding. Although the explanation have been exemplified to the 
Coi8Mn2 cluster, it is applicable to the rest of calculated clusters with 1 < n < 8. For these 
clusters, the DOS of Mn atoms shows a higher-lying occupied bonding orbital for spin-up 
electrons and a lower-lying (a few eV above HOMO level) empty antibonding orbital for 
spin-down electrons [see Fig. [2](d)], while for Co first neighbor atoms the DOS manifests an 
increase of the spin-up population [see Fig. [2]^b)]. In consequence, and based on the Mulliken 
population analysis, the electrons involved in the Co-Mn bonding move from the spin-down 
channel of Mn atoms to the spin-up channel of the Co first neighbor atoms. This mechanism 
contributes to raise the magnetic moment of both Mn and Co neighbor atoms [see inset and 
text of the caption in Fig. ([5])]. However, the process stops at a critical concentration of 
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the Mn atoms because the spin-up channel of the available Co atoms does not admit more 
electrons. Thus, in the range with n varying from 9 to 19, the process reverse and there is a 
charge transferring from the Mn spin- up channel to the Co neighbors spin-down one. In this 
case, the Mn and also the Co neighbor atoms [see inset in Fig. (^] decrease their moments. 

It is also interesting finally to comment that the non-monotonic decrease of the magnetic 
moments at the area above 40 % of Mn concentration (see Fig. [5]) is attributed to the 
discontinuity of g(n) with the concentration of the Mn impurity. Thus, for example, g(n) is 
equal to 1 in the range 9 < n < 12. The average magnetic moments of the AF Mn atoms 
reported in Table IIVI are approximately constant in this range and consequently the total 
magnetic moment per atom of the clusters with 9 < n < 12 is not reduced. The behavior 
changes drastically for n=13, where g(13)=3 and fil^^_ ^^=-2. 39 ^B/(itom. The onset of the 
discontinuity in g(13) and the consequently enhancement of reduction of 

the total magnetic moment of CoyMuis cluster. The same kind of explanation still persists 
for the rest of clusters with higher Mn concentration. 

IV. SUMMARY 

In summary, the magnetic moment enhancement of CoMn clusters below 40% of Mn 
concentration and its successive fast dropping has been explained using first-principles elec- 
tronic structure calculations. The explanation resides mainly on the magnetic role played 
by the Mn atoms and the "spin-flipping" of the electrons belonging to the Co-Mn bonding. 
Moreover, a new formula [Eq. ([5])] has also been proposed for calculating the magnetic mo- 
ments of the CoMn alloy clusters. This equation represents an improvement over the results 
provided by both the RB model and the VBS approximation. 
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where 6{n) is the Heaviside theta function 
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FIG. 2: (Color online) Plot of the spin-polarized DOS for (a) Coi8Mn2, (b) Co-1 and (d) Mn-4, 
where the labels 1 and 4 represent the site number of Co and Mn atoms, respectively. In (b) 
and (d), we have only plotted the contribution of the d orbitals because s and p orbitals are of 
less importance. The upper-half part of the figures is for spin-up electrons while the lower-half 
is for the spin-down electrons. The atoms of (b) and (d) are labeled in the geometric structure. 
We have also plotted the shape of the delocalized molecular orbitals for the higher-lying occupied 
and lower-lying unoccupied levels of Coi8Mn2 cluster with the aim of showing the bonding and 
antibonding character of the orbitals, respectively. The dotted vertical lines represent the HOMO 
level. In (c), DOS at the HOMO level is plotted against the Mn concentration. 
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FIG. 3: (Color online) Illustration of the ground state structures for CoMn binary clusters with 
0, = 20. Co atoms are shown with yellow spheres whereas Mn atoms are represented with spheres 
in magenta. The green color represents the Mn atoms that couple antiferromagnetically with the 
rest of atoms of each cluster. 
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FIG. 4: (Color online) Chemical order parameter as a function of the Mn concentration. The 
solid line is a polynomial fitting to the numerical values of T and serves as a guide to the eye to 
appreciate the asymmetry of the chemical order parameter with respect to the dash vertical line 
(midpoint of the Mn concentration). 
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FIG. 5: (Color online) Magnetic moments per atom of ComMn„ clusters with = 20 calculated in 

led star symbol). The data for the CoMn 



this work (empty star symbol) and measured in Ref. y (fi^ 



bulk are taken from Ref. 



(O symbol). The solid line represents 



21| (A and □ symbols) and Ref. l22 
the linear fitting of the magnetic moments for the clusters with n<8 and it retains a slope of 1.83 ^b- 
The diamond points are the expecting magnetic moments according to Eq. ([SJ . The inset represents 
the difference between the magnetic moments per atom of the Co first neighbors of Mn atoms 
(rhaiioy) in the ComMn„ and their Co counterparts (mco2o) in C020 cluster (Am = rhaiioy ~'^Co2o) 
as a function of Mn concentration. The linear fitting is a guide to the eye. 
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FIG. 6: (Color online) Isosurface plot of the magnetization density (m(r) = [/5^(r) — p^{r)]fiB) for 
(a) CoiiMng and (b) Coi2Mn8. The Co and Mn atoms are illustrated by the yellow and magenta 
spheres, respectively. The blue surfaces represent a positive value for the magnetization density 
whereas the red surface indicates a negative value. 
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FIG. 7: (Color online) Plot of the spin-polarized DOS for (a) Mn-13 and (b) Mn-16 atom of 
CoiiMng cluster. The upper-half area of figure (a) and (b) is for spin-up electrons while the 
lower-half area is for spin-down electrons. The dotted vertical line represents the HOMO level. 



20 



